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Abstract 

We describe the physical hypothesis in which an approximate model of water waves is 
obtained. For an irrotational unidirectional shallow water flow, we derive the Camassa- 
Holm equation by a variational approach in the Lagrangian formalism. 

1 Introduction 

The Camassa-Holm equation reads 

U t + 2kU x + 3UU X - U txx = 2U X U XX + UU XXX (1) 

with x £ R, t € R, U(x,t) € R. Subscripts here, and later, denote partial 
derivatives. The constant k is related to the critical shallow water speed. For 
k=0 the equation (Q]) possesses peaked soliton solutions ([2]). The physical 
derivation of as a model for the evolution of a shallow water layer under 
the influence of gravity, is due to Camassa and Holm [3]. See also Refs. [5], 
|16j for alternative derivations within the shallow water regime. In [19 it was 
shown that the equation (Q} describes a geodesic flow on the one dimensional 
central extension of the group of diffeomorphisms of the circle (for the case 
k—0, a geodesic flow on the diffeomorphism group of the circle, see also [9]). It 
should be mentioned that, prior to Camassa and Holm, Fokas and Fuchssteiner 
[12] obtained formally, by the method of recursion operators in the context of 
hereditary symmetries, families of integrable equations containing ([1]). These 
equations are bi-Hamiltonian generalizations of the KdV equation and possess 
infinitely many conserved quantities in involution (|12j). In order to see how 
the equation ([I]) relates with one of equations in the families introduced in [12] , 
see [14], §2.2. Also, the equation (5.3) in [13] is the equation (JlJ but with errors 
in the coefficients. 
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The Camassa-Holm equation attracted a lot of interest, due to its complete 
integrability [10] (for the periodic case), [8] and the citations therein (for the 
integrability on the line) , the existence of waves of permanent form and of break- 
ing waves [7] and the presence of peakon solutions [3] (for a more rigorous study 
on the peakons see for example [2]). An important question to be answered 
is whether solutions of the water wave problem can really be approximated by 
solutions of the Camassa-Holm equation. In [18] it is shown that suitable solu- 
tions of the water wave problem and solutions of the Camassa-Holm equation 
stay close together for long times, in the case e = S 2 , e being the amplitude 
parameter and 5 the shallowness parameter. In [TH] it is also shown that the 
peakon equation cannot strictly be derived from the Euler equation and hence 
it is at most a phenomenological model. 

The present paper is concerned with the physical hypothesis in which an ap- 
proximate model of water waves is obtained, with the derivation of the Camassa- 
Holm equation by a variational approach in the Lagrangian formalism, and with 
the role of this equation within the shallow water problem. 

2 The governing equations for water waves 

We consider water moving in a domain with a free upper surface at z = 
ho + r/(x,t), for a constant ho > 0, and a flat bottom at z = 0. The undis- 
turbed water surface is z = Jiq. Let (u(x, z,t),v(x, z, t)) be the velocity of the 
water - no motion take place in the y-direction. The fluid is acted on only by 
the acceleration of gravity g, the effects of surface tension are ignored. For the 
gravity water waves, the appropriate equations of motion are Euler's equations 
(EE )(see [15]). Another realistic assumption for gravity water wave problem is 
the incompressibility (constant density p) (see |15j). which implies the equation 
of mass conservation (MC). The boundary conditions for the water wave prob- 
lem are the kinematic boundary conditions as well as the dynamic boundary 
condition. The kinematic boundary conditions (KBC) express the fact that the 
same particles always form the free water surface and that the fluid is assumed 
to be bounded below by a hard horizontal bed 2 = 0. The dynamic boundary 
condition (DBC) express the fact that on the free surface the pressure is equal 
with the constant atmospheric pressure denoted po- Summing up, the exact 
solution for the water-wave problem is given by the system 



u t + uu x + vu z = ~Vx 
v t + uv x + vv z = -- p Vz ~ 9 



(EE) 



Ux + v z = 



(MC) 



(2) 



v = Vt + urj x on z = ho + r](x, t) 
v = on z = 



(KBC) 



p = po, on z = h +T)(x,t) 



(DBC) 
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where p(x, z,t) denotes the pressure. 

We non-dimcnsionalise this set of equations and boundary conditions using 
the undisturbed depth of water, /io, as the vertical scale, a typical wavelength 
A, as the horizontal scale, and a typical amplitude of the surface wave a (for 
more details see [TS], [TB]). An appropriate choice for the scale of the horizontal 
component of the velocity is y/gho. Then, the corresponding time scale is 



fgha 

and the scale for the vertical component of the velocity is ho ^ 9 ^° . Thus, we 
define the set of non-dimensional variables 

± nr, f h + 

(3) 

where, to avoid new notations, we have used the same symbols for the non- 
dimensional variables x, z, 77, t, u, v, in the right-hand side. The partial deriva- 
tives will be replaced by 




(4) 



Let us now define the non-dimensional pressure. If the water would be station- 
ary, that is, u = v = 0, from the first two equations and the last condition with 
77 = 0, of the system ([2]), we get for a non-dimensionalised z, the hydrostatic 
pressure po + pgho(l — z). Thus, the non-dimensional pressure is defined by 

P >-> Po + pgho(l - z) + pgh p (5) 

and 

Px >-> P—Px, Pz >-> -P9 + A9Pz (6) 

Taking into account ([3]), |@]), (JSJ), and ([6]), the water-wave problem ([2]) writes 
in non-dimensional variables, as 

u t + uu x + vu z = —p x 
S 2 (v t + uv x + vv z ) = -p z 

u x + v z = (7) 
v = e(rjt + urj x ) and p = er/ on z = 1 + £77(2;, t) 
v = on z = 



where we have introduced the amplitude parameter e = and the shallowness 

parameter £ = ^ . The small-amplitude shallow water is obtained in the limits 
e — * 0, 6 — > 0. We observe that, on z = 1 + e?7, both i> and p are proportional 
to e. This is consistent with the fact that as e — > we must have t) — ► and 
73 — > 0, and it leads to the following scaling of the non-dimensional variables 



p h-> ep, (it, u) 1 — ► e(it, -u) 



(8) 
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where we avoided again the introduction of a new notation. The problem |(7J) 
becomes 

u t + e(uu x + vu z ) = -p x 
8 2 [v t + e(uv x + vv z )] = -p z 

u x + v z = (9) 
v = rj t + eurj x and p = r) on z = 1 + (rj(x, i) 
v = on z = 

Further, the parameter S can be removed from the system (0) (see [H]), this 
being equivalent to use only ho as the length scale of the problem. In order to 
do this, the non-dimensional variables x, t and v from ([3]) are replaced by 

Je J~e 5 . 

X I ► — —X, t I ^ —=~t, V I — > -p!) (10) 

Therefore the equations in the system ([9]) are recovered, but with S 2 replaced 
by e in the second equation of the system, that is, this equation writes as 

e[v t + e(uv x + vv z )} = -p z (11) 



3 The Camassa-Holm equation 

The non-dimensionalisation and scaling presented above will be useful in ob- 
taining a scheme of approximation of the governing water-wave problem ([2]) . 

In the limit e — > 0, that is for small-amplitude waves, the system ([9]) with 
the second equation given by (jTTJ) becomes 

Ut + Px = 
p z =0 

u x + v z = (12) 
v = rjt and p — r\ on z = 1 
v = on z = 



From the second equation in (fT2")) . we get that p does not depend on z. Because 
p = rj{x, t) on z = 1, we have 

p = rj(x, t) for any < z < 1 (13) 

Therefore, using the first equation in (p~2|) . we obtain 

u = — J r] x (x,t)dt + !F(x,z) (14) 

where is an arbitrary function. Differentiating (|14p with respect to X and 
using the third equation in (|12p . we get, after an integration against z, 



v = z r) xx (x, t)dt - Q(x, z) + g(x, 0) 



(15) 
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where Q z (x, z) = J- X (x, z) and we have also taken into account the last condition 
in the system (fl"2]) . Making z = 1 in (fT5|) . and taking into account that v = r\ t 
on z = 1, we get after a differentiation with respect to t, that 77 has to satisfy 
the equation 

Vtt - Vxx = (16) 

The general solution of this equation is r](x,t) — f(x — t) + g(x + t), where / 
and g are differentiable functions. It is convenient first to restrict ourselves to 
waves which propagate in only one direction, thus, we choose 

r,(x,t)=f(x-t) (17) 

Therefore, for u and v in (fl4]) . (p~5|) we get 

u = X) + T{x,z), v = —zr] x — Q(x, z) + G(x,0) (18) 

with Q z {x,z) = T x {x,z), Q{x, 1) = Q(x, 0), arbitrary functions. Thus, the 
solutions to the shallow water problem are determined by the evolution of the 
function r)(t,x), which represents the displacement of the free surface from the 
undisturbed (flat) state. 

Under the assumption that the fluid is irrotational, we get 

T{x,z) = const := cq and Q(x,z) = (19) 



Indeed, if the fluid is irrotational the vorticity is zero, that is, in addition to the 
system ©, we also have the equation 

u z - v x = (20) 

For a discussion of the role of vorticity in water wave flows see for example [11] , 

El- 

In the equation (|20|) . the velocity components u and v are written in the physical 
(dimensional) variables. If we non-dimensionalise this equation using ([3]), (U]), 
we obtain 

u z = S 2 v x (21) 
After scaling ([5]) and transformation (|10p . the equation (f2"Tj) writes as 

u z = ev x (22) 

Therefore, in the limit e — > 0, we get in addition to the system (fT2")l . the equation 

u z = (23) 

The relation (flU)) remains the same but instead of (|14[) we have now 

u=- r) x (x, t)dt + T(x) (24) 
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where T is an arbitrary function. Using the third equation in (|12[) , we get now 

v = -zu x = z(j r) xx (x,t)dt- P{x)j (25) 

where we have taken into account the last condition in the system (fT2")l . Making 
z = 1 in (|25[) . and taking into account that v = rj t on z = 1, we get after a 
differentiation with respect to t, that ij has to satisfy the equation (TT6")) . We 
consider the solution of (jTHJ) into the form the solution (fT7|) . Therefore, for it 

and w in (pM|) . (|23|) . we have u = ?y + .F(a;), u = — z (rj x + !F'(x)J . The condition 

v = r\t on z = 1, yields .F(a;) = const := en. Thus, for the irrotational case the 
solution of the system (TT2"]) plus the equation (|23p , can be written into the form 

r)(x,t) = f(x -t), u = r) + c , v = -zr) x (26) 

We underline the fact that in our approximation the vertical velocity component 
maintains a dependence on the z-variable. In analyzing the motion of the fluid 
particles, this means that the particles below the surface may perform a vertical 
motion. This is in agreement with a recent general result obtained in \§\ for the 
Stokes waves, i.e. particular waves in an irrotational flow which are solutions of 
the full Euler equations. 

By consistently neglecting the e contribution, we will derive using variational 
methods in the Lagrangian formalism (see [S]), the equation (TTJ governing uni- 
directional propagation of shallow water waves. 

In the Lagrangian formulation of a fluid, the flow pattern is obtained by de- 
scribing the path of each individual water particle. Consider the ambient space 
M whose points are supposed to represent the fluid particles at t = 0. A diffeo- 
morphism of M represents the rearrangement of the particles with respect to 
their initial positions. The motion of the fluid is described by a time-dependent 
family of orientation-preserving diffeomorphisms -f(t, •) £ Diff(Af). A point x in 
M follows the trajectory 7(t, a:) through M. For our problem, since a particle 
on the water's free surface will always stay on the surface and describes pro- 
gressive plane wave (no motion take place in the y direction), we may regard 
the motion at that of a one-dimensional membrane. For the one-dimensional 
periodic motion, M — S 1 the unit circle. We can allow M — R and add the 
technical assumption that the smooth functions defined on R with value in R 
vanish rapidly at ±00 together with as many derivatives as necessary (see [3] 
for a possible choice of weighted Sobolev spaces). In what follows we focus on 
the latter situation. 

For a fluid particle initially located at x, the velocity at time t is 74 (t, a;), 
this being the material velocity used in the Lagrangian description. The spatial 
velocity, used in the Eulerian description, is the flow velocity w(t,X) = j t {t,x) 
at the location X = 7(4, a;) at time t, that is, w(t, •) = 7t 7 _1 - In Lagrangian 
description, the equation of motion is the equation satisfied by a critical point of 
a certain functional (called the action) defined on all paths {j(t, •), t € [0, T]} in 
Diff(R), having fixed endpoints. Following Arnold's approach to Euler equations 
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on diffcomorphism groups (PQ), the action for our problem will be obtained by 
transporting the kinetic energy to all tangent spaces of Diff(R) by means of 
right translations. For small surface elevation, the potential energy is negligible 
compared to the kinetic energy. Taking into account (|26|) , the kinetic energy on 
the surface is 



1 



1 



K = ^ / {u 2 + v 2 )dx = ^ I [u 2 + (1 + er,) 2 u 2 x ] dx « ^ / (u 2 + u 2 x ) dx 

z J -oo L J -oo L J -oo 

(27) 

to the order of our approximation (see [5]). We observe that if we replace the 
path -f(t, •) by -f(t, •) o for a fixed time-independent tp in Diff(R), then 

the spatial velocity is unchanged 74 o 7 _1 . Transforming K to a right-invariant 
Lagrangian, the action on a path j(t, •), t £ [0, T], in Diff(R) is 



0(7) 



1 

2 7o 



de 



{{ito^- 1 ) 2 + [d x { lt o 1 -t)] 2 }dxdt 



0. 



(28) 



The critical points of the action (|28|) in the space of paths with fixed endpoints, 
verify 

d 1 

— a( 7 + e^) 



e=0 



(29) 



for every path ip(t, •), £ e [0, T], in Diff(R) with endpoints at zero, that is, 
<p(0, •) = = <ys(T, •) and such that 7 + eip is a small variation of 7 on Diff(R). 
Taking into account (|28|) . the condition (|29|) becomes 



l 7t 07 



e=0 



[(7t + eVt)° h + s<p) x ] 



+d x (j t o 7 _1 )4 [3*((7t+^)°(7 + £<p) l )]jda;* = 

u£ e=0 J 

(30) 



e=0 



After calculation (see [5]), we get 
d 



de 



e=0 



[{-/t + eipt) o (7 + eip) x ] 



and 



de 



e=0 



[d x (( 7 t + eip) o (7 + e<p) x )] 



^07 1 -(( /3 o7 1 )d x (^ t ° r y x ) 

<3 t (</J o 7 _1 ) + ( 7i o 7 _1 )a i (<^ o 7 _1 ) 
-(^o 7 - 1 )^(7 t o 7 - 1 ) (31) 

d x (<Pt 7 _1 ) - 9a.(7t o 7~ 1 )5 X (</? o 7 _1 ) 
-(^o 7 - 1 )^(7 t 07- 1 ) 

3 te (( / 3 7 - 1 ) + ( 7t 7 - 1 )9^0 7 - 1 ) 

-(^o7- 1 )^(7 t o 7 - 1 ) (32) 
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where are used the formulas of the type 



d_ 

de 



e=0 



d x {lt°l X ) 



w o 7 
lx ° 7" 
l 



Tfa 7 



(33) 



(34) 



7a; ° 7" 

dt(v°7 -1 ) = tp t oj~ 1 + ((p x o 7 _1 ) 3 t (7 _1 ) = ^ f o7 _1 -(7 t o7)3 ;1 .(( y 907~ 1 ) (35) 
Thus, denoting j t ° 7 _1 = u, from (f3"Tj) . ([32]) , the condition (f30|) writes as 

f T /-oo 

/ {u [9 t (y> o 7 _1 ) + ud x (<p o - ((po ■y~ 1 )u x ] 

/ — oo 

+U X [dtx{<P ° 7 _1 ) + U-dliV ° 7 _1 ) - (</3 7 _1 ) M a;a;] } ^(ft = 

(36) 

We integrate by parts with respect to t and x in the above formula, we take 
into account that ip has endpoints at zero, the smooth functions defined on R 
with value in R, together with as many derivatives as necessary, vanish rapidly 
at ±00, and we obtain 



(cp o 7 1 )[u t + 3uu x - u txx - 2u w u a 



ixxx] dxdt = (37) 



Therefore, we get that for an irrotational unidirectional shallow water flow, the 
horizontal velocity component of the water u(x,t) satisfies the Cammassa-Holm 
equation (QJ) for k = 0. 



Let us see now which equation fulfill the displacement r](x, t) of the free 
surface from the flat state. Taking into account (|2"6"|), the kinetic energy on the 
surface is 



2 



(u 2 +v 2 )dx 



[(r/ + c ) 2 + (l + er,) 2 r] 2 x ]dx 



dx 



(38) 



to the order of our approximation. Transforming if to a right-invariant La- 
grangian, the action on a path T(t, •), t € [0, T], in Diff(R) is 



1 f T f°° 

a(r) = 5 / / {(T* o r- 1 + CQ ) 2 + [s a (r t o r- 1 )] 2 }^ 

Z JO J -oo 



(39) 



The critical points of the action (|39[) in the space of paths with fixed endpoints, 
verify 

d -a(T + e^) =0, 



de 



e=0 



(40) 
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for every path <£>(£, •), t € [0, T], in Diff(R) with endpoints at zero, that is, 
$(0, •) = = $(T, •) and such that T + e$ is a small variation of T on Diff(R). 
Taking into account (|3"9")l . the condition PO)) becomes 



T poo 



e=0 



[(r t + e$ t )°(r + e$) -1 ] 



+d x (T t o T- 1 )— [d x ((T t + e$ t ) o (T + e*)- 1 )] ^ dzcft = 

tt£ e=0 J 

(41) 

After calculation, denoting T f o = r), (j^Tj) writes as 

fT /-oo 

{(77 + c ) o r- 1 ) + r?s x ($ o r- 1 ) - (* o r- 1 )^] 
+^ [9tx(# o r- 1 ) + r^($ o r- 1 ) - (* o r- 1 )^,] } dxdt = o 

(42) 

We integrate by parts with respect to t and x in the above formula, we take into 
account that $ has endpoints at zero, the smooth functions defined on R with 
values in R, together with as many derivatives as necessary, vanish rapidly at 
±oo, and we obtain 



($oT x ) [r) f + 'irjrj x + 2c a i] x - r\ txx - 2i] x i] xx - i]i] xxx } dxdt = (43) 



Therefore, we get that for an irrotational unidirectional shallow water flow, the 
displacement r)(x, t) of the free surface from the flat state, satisfies the Camassa- 
Holm equation {!]) for k = cq. 
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